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It was recently pointed out that simple scaling properties of Polyakov correlation functions of 
gauge systems in the confining phase suggest that the ratios of k-string tensions in the low temper- 
ature region is constant up to terms of order T^ . Here we argue that, at least in a three-dimensional 
Z4 gauge model, the above ratios are constant in the whole confining phase. This result is ob- 
tained by combining numerical experiments with known exact results on the mass spectrum of 
an integrable two-dimensional spin model describing the infrared behaviour of the gauge system 
near the deconfining transition. 
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1. Introduction 

It was recently pointed out that in most confining gauge theories, besides the fundamental 
string (of tension a) which is formed between a pair of static sources in the fundamental represen- 
tation /, there is the freedom of taking the sources in any representation If, for instance, the 
gauge group is SU(A^) there are infinitely many irreducible representations at our disposal. How- 
ever, as the sources are pulled apart, no matter what representation is chosen, the asymptotically 
stable string tension depends only on the A^— ality k of M. As a consequence the heavier strings 
decay into the string of smallest string tension Ou- The corresponding string is referred to as a 
k-string. This kind of confining object can be defined whenever the gauge group admits more than 
one non trivial irreducible representation. 

In a previous work [|l]], some of us have proposed an expression for the low temperature asymp- 
totic expansion for these string tensions. An interesting consequence of such an expansion is that 
their ratios are expected to be constant up to terms. The low temperature data presented in sup- 
port of this expectation were taken from Monte Carlo simulations on a particular system, namely a 
(2-i-l)-dimensional Z4 gauge model. 

The main conjecture we want to verify in this work is that OkiT) / o{T), at least in that Z4 
gauge system, is in fact independent of the temperature in the whole of the confining regime. To 
check this idea we used the fact that the Svetitsky-Yaffe (SY) conjecture allows to reformulate 
the system in a totally different perspective, based on a two-dimensional integrable theory. 

It turns out that the deconfinement transition of the 3D Z4 gauge model is second order and, 
according to the SY conjecture, belongs to the same universality class of the 2D symmetric Ashkin- 
Teller (AT) model. 

The two-dimensional AT model can be seen in the continuum limit as a bosonic conformal 
field theory plus a massive perturbation driving the system away from the critical line (i. e. a sine- 
Gordon theory). Thus, a map between the AT critical line and the sine-Gordon phase space is 
provided. This theory is integrable, and the masses of its lightest physical states (first soliton and 
first breather mode, of masses M and M\) correspond to the tensions o{T) and 02{T) near T^, 
whose ratio, in this context, can be analytically evaluated and turns out to be 

lim^ = ^=2sin^(2v-l), (1.1) 
T^T, o{T) M 2^ ' ' ^ ^ 

where v is the thermal exponent. 

1.1 The (2+l)d Z4 gauge model and its dual reformulation 

The most general form of Z4 lattice gauge model admits two independent coupling constants, 
with partition function 

Z(/3/,%) = n I eUPf%+kf%f/^+^-\ (1.2) 

in which the gauge field Ui on the links on a cubic lattice is valued among the fourth roots of the 
identity and the sum in the exponent is taken over the elementary plaquettes of the lattice. Such a 
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theory can be reformulated as two coupled Z2 gauge systems (see for details): 

mf,M = I\ I ,{Up = Y\Ui-Vp = Y\Vi). (1.3) 

/ {c/;=±i,y/=±i} /ep i&p 

From the data in |[T]], obtained by means of finite-temperature measurements of Polyakov- 
Polyakov correlation functions, and particularly from those referring to the point P identified 
by (o;,j8) = (0.050,0.207), the string tensions a and 02 can be evaluated in the T — > limit as 
temperature-independent quantities : 

aa^ = 0.02085(10) , 020^= 0.03356(22) , 

where a is the lattice spacing. Their ratio, which has been argued to equate the central charge of 
the CFT related to the 2-string, is then given by 

— = 1.610(13). (1.4) 



2. The Svetitsky-Yaffe conjecture and the Sine-Gordon model 

The mapping induced by the Svetitsky-Yaffe conjecture leads to a substantial simplification in 
the study of the critical properties of the deconfining transition, allowing to study it as a standard 
symmetry-breaking transition which takes place in a spin model. In the present case we deal with 
the symmetric Ashkin-Teller model in two-dimensions. 

The action for this model is given by: 

^Ar = - £ [J{olal + olo^) +74(aia»,2)] . (2.1) 

Such a model has been extensively studied in the past, and a number of exact results have been 
derived [|3|]. It is useful to note that it can be seen as a perturbation of the Gaussian model, and in 
such a bosonic language the thermal perturbation can be written as cosj3(p, where j8 is a marginal 
parameter equivalent to J/^. Hence we are left with 

jT^^ = y Q5^(p(9^(p - Tcosj8<pj , (2.2) 

which is the action of the Sine-Gordon model. Furthermore, since the confined phase of the gauge 
theory is mapped in the high-T phase of the Ashkin-Teller, we will only consider the case T > 0. 

Such a QFT is of particular interest because it is integrable, and this is the main reason for 
rewriting the action of the model near the critical point in a bosonic form. Integrability means that 
an infinite number of integrals of motion exists. The main consequence in (l-i-l) dimensions is the 
fact that the scattering theory is very constrained, because the S-matrix is factorised in products 
of two-body interactions, and inelastic processes are forbidden. It follows that the latter can be 
computed exactly together with the mass spectrum. 
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2.1 Operator correspondence, mass spectrum and correlation functions 

We already know, from the Svetitsky-Yaffe original work, that the Polyakov loop in the funda- 
mental representation corresponds to the spin operator. Then, following the same reasoning used 
in [Q], it is possible to deduce that the Polyakov loop in the double fundamental representation is 
related to the so-called polarisation operator ^ = a'a^, where a' and are the spin variables 



defined in (2. 1 ). Its bosonic form and the corresponding anomalous dimensions are given by 



^ = sin|(p, ^^ = ^; (2.3) 

we also notice that (^) = in the high-T phase of the model. 

Sine-Gordon mass spectrum [|[] : The exact knowledge of the S-matrix allows to access to the 
exact mass spectrum of the theory. Without entering the details, the spectrum of the SG model 
is given by a soliton/anti-soliton doublet of fundamental particles of mass M, and a number of 
soliton/anti-soliton bound states, called breathers B„, whose number is a function of jS^. By defin- 
ing the coupling constant ^ in the following way 

E = ^ (24) 
^ 871-/32' 

we have that for ^ > n, i.e. > 4n, no bound states are present and hence the spectrum is given 
by the soliton/anti-soliton doublet only (repulsive regime). For ^ < 7i, i.e. < 4n, we are in the 
attractive regime and the breathers B„ appear as simple poles of the S-matrix (see for example 
The next step is to associate particle states to operators in the high temperature phase. It has been 
done in by taking into account their properties of symmetry and locality. The consequence is 
that the spin operator is naturally associated to the mass of the soliton, and the polarisation operator 
is associated to the mass M\ of the breather Si. Hence, following the Svetitsky-Yaffe conjecture, 
the ratio of string tensions in the confining phase near the transition is given by 

-=2sm-. (2.5) 

This result, being a dimensionless ratio, is expected to be universal in the limit T ^ 0. This fact can 
be explicitly seen by expressing the coupling ^ in terms of some critical exponent. It is possible to 
work out the following relation between ^ and the thermal critical exponent v 

^=n{2v-l) ^ ^ = 2sin|(2v-l). (2.6) 

Correlators at large distance: The previous analysis of the mass spectrum allows to compute 
the leading behaviour of the correlators (era) and {^^) at large distance by means of their spec- 
tral expansion over form factors (the interested reader can refer to |8|] for the details). 

The analysis of the previous section allows immediately to write down the leading term for 
(aa) and (^^) correlators in the high-T phase of the theory, up to an inessential proportionality 
constant 

{a{x)a{0)) ^ Ko{M\x\), \x\^^; 
{^ix)^{0)) KoiMi\x\), W^oo, (2.7) 

where ^0 denotes the modified Bessel function of order zero, and M, Mi are the masses of the 
sohton and the first breather respectively. 
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2.2 Baryon vertices and mass spectrum 

As noticed in [|l|], the balance of the string tensions for a given vertex gives the following 
expression for the angles at the center of the junction of three arbitrary k-strings 

a^{T) + al{T)-of{T) 

cos di = — — - — - — , and cyclic permutations of the indices. (2.8) 

20j{T)0k{T) 

The rigidity of the geometry of the vertex is then ensured by requiring that such angles are kept 
fixed when the temperature varies. As a consequence, all the string tension ratios are constant up 
to a given order in T , namely as far as the effective string picture is valid. 

A similar picture emerges when studying the gauge system near the deconfining transition. 
The scattering theory describing the system in such a case can exhibit bound states whose mass 
is given by the following relation 

= + 2/721^2 cos Mj2, triangle of masses (2.9) 

where Q = / purely imaginary value of the rapidity corresponding to the creation of the 

particle mi,, and m\, m2 are the masses of the initial state. 

In the present case the process of coalescence of two fundamental strings into a 2-string corre- 
sponds to the scattering of a soliton/anti-soliton pair creating the bound state B\. For such a process 

Bi 

^ss 



we know that m^^ = k — £, which, once inserted in (2.9), gives 



M^ = 2M2(1-coS(§) ^ ^=2sin| (2.10) 



which is nothing but the mass formula used in the previous Section. 

3. Monte Carlo setting and procedure 
3.1 Mass ratio by correlators 



As introduced in Subsection we can determine the ratio M\ /M using the large distance 
asymptotic behaviour of correlators; actually, exploiting the Svetitsky-Yaffe conjecture, we mea- 
sured the Polyakov-Polyakov correlators Gag{R) of the (2-i-l)d Z4 gauge theory: 

G,*W = (P«(0)PiW). (3.1) 

In Section |TTT] we have explained we can study this theory by means of simulations on the AT 
model and in the measurement of Polyakov-Polyakov correlators in both the fundamental and 
double fundamental representations, G{R)f and G{R)ff, is described in detail. 

We have taken 10^ measures on the 64^ x 7 lattice in the phase space point P;Nx = l is chosen 
because it is the lowest possible value above the deconfinement transition. Simulations have been 
done for each value of R in the range [15 -^44]. These data are fitted using an expansion of the 
Ko{mR) Bessel function, truncated to first two terms, in a range [Rmin,Rmax\, where Rmax = 44; we 
have verified the results are stable when /?,„„, varies in the range [22 33] . Therefore, it is possible 
to determine the two masses: 

aMff = 0.0698(15) {x^/d.o.t w 1.3) , aMf = 0.0433(8) (x'^/d.o.f. w 1.2) , 
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from which we can determine the ratio: 

Oi{T ~ Tc)/o{T ~ T,) = Mff/Mf = 1.612(46) . (3.2) 
This result, obtained near the critical temperature, is compatible with the zero-temperature value 



(1.4), providing a strong evidence for our conjecture. 



3.2 Estimating 02/0 through the thermal exponent v with finite-size scaling 



To use the formula for the mass ratio, Eq. (2.6), we need a quite precise estimate for the thermal 



critical exponent v in the phase space point P. It can be obtained by means of a finite-size scaling 
analysis of the plaquette operator or some related observable that we denote with ( □ )l, where L is 
the spacial size of the lattice. Let us notice that in the SY context the plaquette operator is mapped 



into a combination of the unity and the energy operator of the corresponding CFT []10|]. 

In order to exploit the computational advantages of the dual transcription of the gauge model, 
it is convenient to evaluate directly the internal energy of the 3D AT model defined by 

(□)L^-3^(5Ar) . (3.3) 
We decided however to use the corresponding (density of) susceptivity 

{x)L^{{a-{n),f),, (3.4) 

whose power-law to compare with has the form 

{x)L = b'-L^-' , (3.5) 

with the advantage that no constant additive terms are present, which could largely spoil the stabil- 
ity of the numerical results. 

At the practical level, the system at the coupling P turns out to be critical for a temperature 
such that 6 < jr <7, hence, having to work with integer inverse temperatures, it is not possible to 
avoid some approximate method. In particular we decided to define two new points, Pj and P(„ at 
which the system is critical for temperatures T = l/l and T = 1/6 respectively, and then, with a 
linear interpolation, construct the corresponding quantity for the original point of phase transition 
P. To perform the simulations, we used a cluster-based nonlocal update algorithm, an adaptation 
of the Swendsen-Wang prescription, which is described in more detail in [0]. 

We used L = 200 finite-temperature lattices to find the couplings corresponding to P^ and Pj, 
and at such critical points we took ^(10^) measurements of the plaquette at 26 values of spatial 
side L, ranging from L = 10 to L = 165. The data fitted very well to the expectation from L = 70 
already, so we could extract two values of the critical index v: 

Vr=i/6 = 0.8004(19) [22] , Vr=i/7 = 0.7942(18) [38] , 

in which the first uncertainty refers to the statistical fluctuations while the second is an estimate of 
the systematic error in the measurement. 



6 



On the ratio of string tensions in the 3D Z4 lattice gauge theory 



P. Grinza 



By linear interpolation along the couplings, the value of v and the (coupling-dependent) criti- 
cal temperature Tc was calculated for the very point P. We found 



Tc 



^ 1.0393(12), v(/') =0.7984(19) [27]. (3.6) 



By plugging it into the formula for the mass ratio (2.6), we obtain the following result for the 
mass ratio: 

Ml 

^(P) = 1.6124(71)[102] , (3.7) 

which is well compatible with the less accurate estimate coming from the quantities in [|l|] and thus 
well supports our conjecture. 

4. Conclusions 

In this paper we studied the ratio of the string tensions 02{T) /o{T) near the deconfining point 
Tc of a 3D Z4 gauge model and compared the result with a general formula which is expected to be 
true near 7 = for a generic gauge theory in three or four dimensions. In this particular case we 
have combined numerical experiments with known exact results of an integrable 2D quantum field 
theory that belongs, according to the Svetitsky-Yaffe conjecture, to the same universality class of 
the critical gauge system. 

An interesting property of the integrable model is that the mass ratio of the two physical states 
of the theory, which should equate the string tension ratio near T^, can be expressed as a simple 
function of the thermal exponent v (see Eq.(^)). We used two different methods to evaluate 
such a ratio, and both the estimates give compatible results which nicely agree with the ratio 02/0 



evaluated at T = (see Eq.(|1.4|)). We then conclude that, at least in this model, the k-string tension 
ratios do not depend on T . 
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